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ABSTRACT 

We report on a Chandra observation of the massive, medium redshift (z = 0.1386) cooling flow cluster Abell 
1068. We detect a clear temperature gradient in the X-ray emitting gas from kT ~ 5 keV in the outer part of 
the cluster down to roughly 2 keV in the core, and a striking increase in the metallicity of the gas toward the 
cluster center. The total spectrum from the cluster can be fit by a cooling flow model with a total mass deposition 
rate of ~ 150 M Q yr" 1 . Within the core (r < 30 kpc), the mass depositon rate of ~ 40 M Q yr" 1 is comparable to 
estimates for the star formation rate from optical data. We find an apparent correlation between the cD galaxy's 
optical isophotes and enhanced metallicity isocontours in the central ^100 kpc of the cluster. We show that the 
approximate doubling of the metallicity associated with the cD can be plausibly explained by supernova explosions 
associated with the cD's ambient stellar population and the recent starburst. Finally, we calculate the amount of 
heating due to thermal conduction and show that this process is unlikely to offset cooling in Abell 1068. 

Subject headings: cooling flows — galaxies: clusters: general — galaxies: elliptical and lenticular, cD — 
intergalactic medium — X-rays: galaxies 



1. introduction 

Galaxy clusters frequently have high X-ray surface bright- 
ness cores due to thermal emission from dense gas. Absent a 
substantial and persistent heat source, this high density, rela- 
tively low temperature gas should cool on a timescale that is 
much less than the age of clusters, leading to a so-called "cool- 
ing flow" (Fabian 1994). The cooling rates reported from the 
low resolution X-ray observatories that operated throughout the 
previous two decades were often extraordinarily large, and fre- 
quently exceeded several hundred to over two thousand solar 
masses per year. Such large cooling rates posed a dilemma. 
They imply the presence of enormous sinks of cold gaseous and 
stellar matter in the cD galaxies found at the centers of cooling 
flows. Although substantial amounts of cold gas and vigorous 
star formation are now found commonly in cD galaxies cen- 
tered in cooling flows (Edge 2001, Edge et al. 2002, Jaffe & 
Bremer 1997, Jaffe, Bremer, & Van der Werf 2001, Donahue 
et al. 2000, Falcke et al. 1998, Crawford et al. 1999, Cardiel, 
Gorgas, & Aragon-Salamanca 1998, McNamara & O'Connell 
1989), they are found in quantities that are much smaller than 
these earlier cooling rates would imply (cf., McNamara 2002). 

This large discrepancy between the cooling and star forma- 
tion rates is one of the primary questions behind what has be- 
come known as the "cooling flow problem". Simply stated, it 
is either a violation of mass continuity, assuming the matter 
is in fact not cooling at the rates implied by the X-ray data, 
or a missing light problem, if the matter is indeed cooling at 
these historically large rates but is hidden in some exotic state. 
Recent high spatial and spectral resolution cluster observations 
made with the Chandra and XMM-Newton X-ray observato- 
ries have changed this canonical dilemma and suggested a third 
possibility. The absence of detected emission lines from gas 
at temperatures less than about 3 keV in XMM-Newton RGS 



and Chandra HETG observations of cluster cores implies that 
a large fraction of the cooling gas does not in fact cool below 
X-ray emitting temperatures (Peterson et al. 2001, Peterson et 
al. 2003, Wise et al. 2004). Instead, this gas is assumed to 
be reheated to ambient temperatures by one or more of several 
possible agents that are now being extensively studied. These 
agents include mechanical and cosmic ray reheating from su- 
pernova explosions, heat conduction from the hot outer layers 
of clusters, and mechanical heating by the central radio sources 
in cD galaxies. 

Images of the cooling flow regions of clusters now show a 
great deal of structure, including the large cavities or bubbles in 
the X-ray emission associated with the radio sources harbored 
by the central cluster galaxies (e.g., McNamara et al. 2000, 
Fabian et al. 2000). The strong interactions between the radio 
sources generated in the nuclei of cD galaxies and the surround- 
ing X-ray emitting gas may reduce the cooling rates and bring 
them more in line with the star formation rates (David et al. 
2001, Nulsen et al. 2002). In addition, Chandra's superb spa- 
tial resolution provides the capability to map the temperature 
structure of the keV gas on fine spatial scales, which permits 
a direct comparison between the sites of star formation seen at 
optical and ultraviolet wavelengths, and sites of rapid cooling. 
In essence, we can now reevaluate the proposition that cool- 
ing flows are fueling star formation by comparing the cooling 
and star formation rates on the same spatial scales. In addition, 
the newly revealed complexity of the gas in cluster cores has 
sparked a revaluation of several possible heating mechanisms, 
such as AGN heating (Brighenti & Mathews 2002, Ruszkowski 
& Begelman 2002), heat conduction from the hotter outer layers 
of clusters (Narayan & Medvedev 2001, Zakamska & Narayan 
2003), and supernovae. 

In this paper, we present a detailed examination of the X- 
ray structure of the core of the Abell 1068 galaxy, and present 
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some of the most detailed temperature, and metallicity maps 
available for a galaxy cluster. In a companion paper (McNa- 
mara, Wise, & Murray 2004, hereafter "Paper 2"), we discuss 
the remarkable starburst properties of the Abell 1068 cD galaxy 
and present a detailed comparison between the X-ray and op- 
tical properties of the cD. In these papers, we show that 90% 
of the ultraviolet light from the starburst is emerging from the 
region of the cluster with the shortest cooling time, and that the 
star formation rates and local cooling rates agree. However, we 
stop short of presenting this fact as proof that the starburst is 
being fueled by the cooling flow, as there are independent data 
to suggest that the starburst may be associated with an ongoing 
dynamical interaction between the cD and a group of compan- 
ion galaxies. In the proceeding analysis, we have adopted a red- 
shift for Abell 1068 of z = 0.1386 (Crawford et al. 1999) and a 
cosmology with Ho=70 km s" 1 Mpc" 1 , n M =03, and S1a=0.7. 
These assumptions yield a luminosity distance of D = 645 Mpc, 
an angular diameter distance of 505 Mpc, and a linear scale of 
2.45 kpc per arcsec. 

2. OBSERVATIONS AND DATA REDUCTION 

Abell 1068 was observed by Chandra on Feb. 4, 2001 for 
30 ksec using the ACIS-S3 back-illuminated CCD (OBSID 
01652). After standard filtering, the net exposure time was 26.8 
ksec. The data were examined for evidence of bad aspect cor- 
rection and background flares but were found to be free of both 
defects. Standard level 2 event files were created using the lat- 
est ACIS calibration products appropriate for the observed focal 
plane temperature of -120 C. All data preparation and creation 
of analysis products (exposure maps, response matrices, etc.) 
employed the CIAO 2.2. 1 data analysis tools and version 2.9 of 
the Chandra calibration database. 

The X-ray morphology of Abell 1068 is complex exhibit- 
ing a highly elliptical appearance on larger scales (~ 300^1-00 
kpc) and complicated structure in the core (r <J 50 kpc). This 
morphology is illustrated in Figure 1 which shows an exposure 
corrected, adaptively smoothed flux image of Abell 1068 in the 
0.3-7.0 keV band. The X-ray emission exhibits a strong cen- 
tral peak centered on the central galaxy characteristic of "cool- 
ing flow" clusters. The exact location of the X-ray peak was 
determined by fitting a two dimensional, elliptical Lorentzian 
surface with varying minor over major axis ratio and position 
angle. The resulting centroid position was (10 40 44.5, +39 
57 12.3) which is within 1.95 arcsec of the optical position of 
the central cD galaxy (10 40 44.62 +39 57 10.2). The best fit 
values for the position angle and axis ratio were 9 = 133 ±3 
and e = 0.71 ±0.02, respectively. These values were used in all 
subsequent spatial and spectral analyses. 

A number of bright point sources are evident in the field of 
Abell 1068 including a bright source ~ 23 arcsec south asso- 
ciated with the radio source B3 1037+401. Although these 
sources contribute a very small fraction of the detected events in 
the field (typically ~ 100-400 counts each), we have excluded 
them prior to extracting the spectral and surface brightness in- 
formation presented here. The CIAO tool WAVDETECT was 
used to determine positions and sizes for point sources in the 
field. These sources were then excised from the event file for 
all further analysis. 

All spectral analysis was done using the ISIS 
(Houck & DeNicola 2000) spectral fitting package which in- 
cludes the standard XSPEC (Arnaud 1996) model library 
as well as access to the high resolution APED database 



(Smith et al. 2001). ACIS effective area files (ARFs) and re- 
sponse matrices (RMFs) were computed using the standard 
CIAO tools and calibration files. A correction for the time- 
dependent degradation of the ACIS effective area due to con- 
tamination buildup on the optical blocking filter was also in- 
cluded. All analysis was conducted using the gain-corrected, 
PI data. A single RMF was computed appropriate for the cen- 
troid position of Abell 1068. Variations in the resolution of 
the CCD, as captured in the response matrices, over the region 
encompassing Abell 1068 were found to be small (<^10%). 
Fits using RMFs computed at other positions in the cluster as 
well as counts-weighted, average RMFs gave equivalent results 
within the errors. Effective area files were calculated at individ- 
ual positions as appropriate. Finally, background files for the 
spectral analysis were extracted from the edge of the ACIS-S3 
CCD, far outside the largest extracted cluster region. 

3. SPECTRAL ANALYSIS 

In order to study the spectral character of the X-ray emis- 
sion from Abell 1068, we have performed three distinct analy- 
ses. First, an integrated spectrum comprising the majority of the 
emission from the cluster was extracted to examine global prop- 
erties such as mean temperature, abundance, and total X-ray 
cooling rate if appropriate. We have also performed a spatially 
resolved analysis using spectra extracted in elliptical annuli to 
study global trends in the temperature, abundance, and mass 
deposition profile as a function of radius in the cluster. Finally, 
we have performed a fully two-dimensional mapping of various 
spectral parameters in the core of Abell 1068 to compare the 
complicated X-ray spatial structure with observed structures in 
other wavebands. We discuss each of these analyses and their 
results here in turn. 

3.1. Integrated Spectrum 

Extracting a "total" spectrum for Abell 1068 is complicated 
by the fact that the cluster spans more than one ACIS CCD. 
While variations in the resolution from different locations on 
the chip are typically small for the S3 aimpoint CCD (where 
the majority of the emission from Abell 1068 is located), they 
can be quite sizable (factors of depending on energy) on 
front-illuminated (FI) CCDs. This difference makes combining 
spectra from ACIS-S3 with spectra from neighboring FI chips 
impractical. Consequently, we have restricted our integrated 
spectral analysis to only that emission contained on ACIS-S3. 
A 175 x 129 arcsec elliptical extraction region centered on the 
X-ray peak and oriented with the measured position angle was 
used (see §2). This region corresponds to a linear extent of 
0.576 Mpc (as measured along the semi-major axis) and com- 
prises ~77% of the total cluster flux as determined from the 
surface brightness analysis discussed below in §4. The resulting 
spectrum contained 65634 counts after background subtraction 
and was adaptively binned to yield a minimum of 30 counts per 
bin prior to analysis. 

We have considered five types of models in fitting the in- 
tegrated spectrum from Abell 1068. In the simplest of these, 
the X-ray emission was represented by a single isothermal 
plasma (hereafter denoted IT) using the MEKAL model from 
the XSPEC library in ISIS where the temperature and elemen- 
tal abundance were allowed to vary. As an intermediate case 
between a single temperature plasma and a fully multi-phase 
model, we also considered a model with two temperature com- 
ponents (2T subsequently). In addition, three models were em- 
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ployed consisting of a single temperature MEKAL model plus a 
multi-phase cooling flow component (MKCFLOW in ISIS). For 
all models featuring a cooling flow, the ambient temperature of 
the MEKAL component and the upper bound on the cooling 
gas were tied. In models with multiple components, the abun- 
dances were tied together and allowed to vary. A foreground 
Galactic absorption column was included in all five models and 
each model was fit twice: once allowing the foreground Galac- 
tic absorption column to vary, and then a second time fixing the 
column to the canonical value of 1 .40 x 10 20 cm" 2 . 

The lower bound on the gas temperature in the cooling flow 
models was chosen to reflect differing scenarios for the fate of 
the cooling material. A scenario where the gas was assumed to 
fully cool below X-ray emitting temperatures was modeled by 
fixing the lower temperature bound at 0. 1 ke V. We have denoted 
this model as FC (full cooling). However, recent XMM RGS 
and Chandra HETG observations of cluster cores have failed to 
detect emission lines from gas at temperatures less than about 
3 keV indicating that cooling may be truncated somehow (Pe- 
terson et al. 2001, Peterson et al. 2003, Wise et al. 2004). This 
situation was replicated by fixing the lower temperature bound 
of the cooling flow component to a value of 2 keV, consistent 
with the observed radial temperature profile discussed in §3.2, 
and has been denoted TC (truncated cooling). Finally, a model 
where the lower limit on the temperature of the cooling gas was 
allowed to be a free parameter was considered (hereafter de- 
noted VC for variable cooling). 

The resulting spectral parameters for all 10 models are pre- 
sented in Table 1 . As a general comment, all models do a rea- 
sonable job at fitting the integrated spectrum, especially at en- 
ergies above 1 .0 ke V, yielding reduced \ 2 values ranging from 
1.1-1.3. Differences in the fits are driven by the spectrum be- 
low 1.0 keV. Models with freely varying absorption do a better 
job than models with the canonical value. Derived columns are 
typically 70% higher than the reported value. This difference 
could be related to ACIS calibration uncertainties below 1.0 
keV. Alternatively, it could represent the presence of excess ab- 
sorption associated with the cooling material or other cold gas 
in the core though this interpretation seems unlikely (see §3.2). 
The single temperature component (IT) fits indicate a global, 
mean gas temperature of ^4 keV and a metallicity of 0.5 so- 
lar abundance. The spectral data, however, clearly indicate the 
presence of multiple temperature components. This result is 
not surprising given the clear temperature gradient in the spa- 
tially resolved analysis. Formally, the best fit is achieved by the 
2T model with temperature components of T\ ^7.0 keV and 
T2 ~ 2.0 keV and a reduced \ 2 = 11. These derived tempera- 
tures are not particularly sensitive to the value of the absorbing 
column. However, this fit results in systematic residuals below 
0.7 keV. The cooling flow models give equivalent values of \ 2 
and better residuals at low energies. 

With CCD resolution spectra alone, it can be difficult to dis- 
tinguish between various cooling flow models. For Abell 1068, 
we find that the full cooling (FC) model with free Galactic ab- 
sorption produces the lowest \ 2 value and most uniformly dis- 
tributed residuals below 1.0 keV. Figure 2 shows the best fit FC 
model overlaid on the data. The resulting cooling rate for the 
FC model was determined to be M x ~ 1 15 — 150 M Q yr" 1 with 
an ambient gas temperature consistent with the global mean. 
This rate drops to ~ 40 M Q yr" 1 if the standard Galactic col- 
umn is used. If the cooling is truncated at 2.0 keV, then the 
resulting deposition rates are increased dramatically reaching 



values of M x ~ 900 — 1000 M Q yr" 1 . Such an increase is to be 
expected of course, since if the gas cools over a smaller tem- 
perature range, then more gas must cool in order to produce 
the same flux as in the FC model. Interestingly, the variable 
cooling models (VC) converge to a preferred lower tempera- 
ture limit for the cooling gas of ~ 1.5 keV with an accom- 
panying small decrease in the cooling rate relative to the TC 
models. This value is consistent with the central temperature 
obtained from our deprojection analysis described below (see 
§3.3. However, the additional free parameters in this model 
render these fits somewhat unstable and make determining con- 
fidence limits difficult. So, while suggestive, one should view 
this result cautiously. 

3.2. Annular Spectral Analysis 

We have examined spatial variations in the physical state of 
the X-ray emitting plasma in Abell 1068 by dividing the ellip- 
tical region containing the previously discussed total spectrum 
into a series of concentric, elliptical annuli. The annuli were 
constructed to have the same axis ratio and orientation angle as 
before and were adaptively sized to include a minimum of 5000 
counts after background subtraction. This requirement yielded 
14 annuli out to a maximum, major-axis radius of 125 arcsec. 
The spectrum outside this radius is dominated by background 
emission. Spectra were extracted from each annular region and 
grouped in the same manner as the integrated spectrum. In 
all subsequent discussion, quoted radii shall refer to distances 
along the cluster's major-axis except where specifically noted. 

In fitting the spectrum from each annuli, we considered the 
various models used to describe the total, integrated spectrum 
in Table 1 . However, the lower signal to noise in the individual 
annular spectra make distinguishing between the various mod- 
els difficult. As a general result, each of the annuli is well fit by 
a single temperate MEKAL plasma (the IT model of Table 1) 
with x 2 values ranging from 0.9-1.3. Additional temperature 
or cooling flow components do not produce significantly bet- 
ter fits. The exceptions to this rule are the central three annuli 
which show a slight improvement in fit quality when a second 
temperature component is included (2T models). Fits for these 
three annuli using the 2T models are consistent with cool gas 
at T ~ 2.0 keV plus a contribution from hotter gas close to the 
global mean value of 4-5 keV. As we discuss in §3.3, the pres- 
ence of these two components is consistent with projection ef- 
fects due to hotter gas along the line of sight. Such contami- 
nation is more easily resolved in the central annuli due to the 
larger temperature difference between the intrinsic cool, core 
gas and the outer, hotter cluster gas. 

There is a clear temperature gradient from 5.5 keV in the 
outer parts of the cluster down to <~ 2.5 ke V in the core (see Fig- 
ure 3). The temperature drops steadily in to a radius of ~ 20 kpc 
(~ 8 arcsec) and then flattens although only two bins cover this 
region. The innermost annulus shows a slight temperature in- 
crease which may be due to emission from either a central AGN 
or the observed starburst in Abell 1068. Fits to this central an- 
nulus including a power-law component to represent emission 
from a central AGN do not however result in a better fit nor 
do they constrain the hardness of the suspected AGN emission. 
For an assumed photon index of 1, we can place an upper limit 
of 5.3 x 10~ 13 ergs s" 1 cm" 2 on the unabsorbed, 0.3-10.0 keV 
flux from any central AGN. This limit represents a maximum 
of ~ 40% of the flux from within the central ~ 4 arcsecs. If a 
photon index of 2 is assumed, this limits rises to 45%. Given 
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these upper limits, the observed temperature increase in the in- 
nermost annulus is more likely to be associated with energy 
injection from the starburst. 

The abundance in the gas shows a steady increase toward the 
center of the cluster from a value of ~ 0.2 solar at the outer 
edge to a central value of ~ 1.0 solar (see Figure 4). There is 
no evidence for a radial dependence on the amount of absorb- 
ing column in the cluster. At all radii, the measured column 
is consistent with the values obtained from fits to the global 
spectrum with free absorption (see Figure 5). This result would 
argue against the presence of intrinsic, excess absorption asso- 
ciated with accumulated cooled gas from a cooling flow. Even 
if ACIS calibration uncertainties make the absolute magnitude 
of the fitted absorbing column suspect, one would expect a rel- 
ative increase toward the center of the cluster if excess absorp- 
tion associated with the cooling gas were present. 

Finally, we have measured the mass deposition rate assum- 
ing a cooling flow model. In contrast to the annular analysis, 
the cooling flow models were fit to concentric elliptical ex- 
traction regions of increasing radius. Therefore, the measured 
value of M x (< r) corresponds to the total cooling rate within 
a given radius. Both the TC and FC models were fit and show 
similar a radial dependence qualitatively. Of the two cooling 
models, the FC model gives marginally better fits although, 
again, the quality of the data does not allow a definitive se- 
lection between the two options. The lower signal to noise of 
the resolved spectra precluded fitting the VC model. Figure 6 
shows the resulting profile for the FC model. In the core region 
(r <; 25 kpc), we find a value of M x (< r) ~ 40 M Q yr" 1 . As 
we discuss in Paper II, this value is consistent with the rates of 
massive star formation implied by optical data in Abell 1068. 
Mx increases steadily with radius reaching a maximum value 
of ~ 150 M yr" 1 in the outer parts of the cluster. Fitting a 
functional form to the mass deposition profile of M x (< r) oc r a 
yields a value of a ~ 0.4, somewhat flatter than the M x (< r)ccr 
relation often seen in pre-Oia/idra cluster data. 

3.3. Deprojected Temperature Profiles 

Due to projection effects, the temperatures measured in 
the previous annular analysis actually represent the emission- 
weighted mean temperature of all gas along the line of sight 
within a given annulus. The spectrum extracted from a given 
elliptical annulus consists of a superposition of emission from 
the local volume bounded by the inner and outer ellipses and 
contributions from hotter gas at all larger radii. Consequently, 
in a standard analysis such as in §3.2, the measured gas temper- 
ature at a given radius will be higher than its actual local value. 
To correct for the effects of this projection, we have performed 
a spectral analysis on the same annular spectra used in the pre- 
vious section accounting for the contribution of hotter gas at 
larger radii to the spectrum from any given annulus. 

To calculate the contributions from gas along the line of 
sight, we have assumed the cluster to be composed of ellipti- 
cal shells with boundaries given by the inner and outer radii of 
the annuli and an axis ratio fixed to the measured value. Since 
we have no information to constrain the inclination angle of the 
cluster, we have assumed the orientation to be perpendicular to 
the line of sight although the equations for the volume contri- 
butions are equivalent for the edge-on case. The intrinsic spec- 
trum of each volume shell was modeled as a single MEKAL 
plasma with variable normalization, temperature, and elemen- 
tal abundance and the fractional contribution of each volume 



shell to the various annuli was then calculated. A foreground 
Galactic absorption component was included and allowed to 
vary. The spectra from all 14 annuli were then fit simultane- 
ously to determine the intrinsic temperature profile in the clus- 
ter. This entire procedure has been implemented as a custom 
model within ISIS (Houck & DeNicola 2000) and is described 
in detail in Wise & Houck (2004). Similar techniques have been 
used by other groups to correct for projection in clusters such 
as Abell 2390 (Allen, Ettori, & Fabian 2001) and Abell 2199 
(Johnstone et al. 2002). 

Qualitatively, the deprojected temperature profile is quite 
similar to that obtained from the standard annular analysis (see 
Figure 7), although with somewhat larger errors. At large radii, 
the temperature agrees with the simple analysis as one would 
expect since contamination from other shells drops with in- 
creasing radius. The biggest change occurs for the central tem- 
perature where we obtain obtain a somewhat cooler minimum 
temperature of 1.8 ± 0.2 keV as compared to the 2.5 keV ob- 
tained in the analysis of §3.2. As noted in the simple annular 
analysis, the innermost annulus remains somewhat hotter than 
the minimum temperature possibly due to contamination from 
an unresolved, central AGN. Values for the abundance were 
consistent with the simple annular analysis and the Galactic col- 
umn was determined to be 2.9 x 10 20 cm~ 2 . 

3.4. Temperature Map 

It is apparent even by visual inspection that the X-ray emis- 
sion from the core of Abell 1068 exhibits a complicated mor- 
phology (see Figure 1). As we discuss in more detail in Paper II, 
the optical data show similar complicated structures including 
condensed knots of extremely blue light presumably associated 
with sites of ongoing star formation (McNamara, Wise, & Mur- 
ray 2004). In order to examine possible correlations between 
the plasma properties and features seen in other wavebands, we 
have constructed 2D maps of the gas temperature and abun- 
dance in Abell 1068. These maps were computed using a cus- 
tom module within the ISIS (Houck & DeNicola 2000) spec- 
tral analysis package which automates the process of extracting 
spectra, generating appropriate ARFs and RMFs, and then fit- 
ting a given spectral model (see (Houck et al. 2004) for a de- 
tailed description). 

The maps were generated using a grid of adaptively sized ex- 
traction regions selected to contain a minimum of 1000 counts 
in the 0.5-7.0 keV band. The grid was selected to span a 
300 x 300 arcsec region centered on Abell 1068. Due to the 
drop in surface brightness with increasing radius, the size of the 
extraction regions increases with radius. Consequently, sharp 
spatial features are not well resolved in the outer parts of the 
maps. In the core (r <; 40 kpc), however, the extraction regions 
range in size from ~ 2-10 arcsecs and features in the temper- 
ature or abundance distributions on these scales can be identi- 
fied. Within each extraction region, the resulting spectrum was 
fit with a MEKAL plasma model including foreground Galactic 
absorption fixed at the nominal value of 2.41 x 10 21 cm~ 2 as de- 
termined from the integrated fits in §3.1. The plasma tempera- 
ture and abundance were allowed to vary and single-parameter, 
90% confidence limits were computed for both parameters at 
each map point. Typical errors for the measured temperatures 
and abundances range from 10% — 30% over the map with the 
most accurate measurements corresponding to the central-most 
regions where the surface brightness is highest. Fits allowing 
the Galactic column to vary produced similar maps within the 
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measured uncertainties and in any event were morphologically 
equivalent. 

// 

The temperature map for the central 80 x 80 of Abell 1068 
is shown in Figure 8. The 0.3-7.0 keV X-ray surface bright- 
ness contours are overlaid for comparison. Although the overall 
temperature structure is consistent with the radial temperature 
profile presented in Figure 3, significant asymmetry in the gas 
temperature is apparent. To the southeast of the cD position, an 
arc of 2-3 keV is present. This extension of cooler gas is spa- 
tially coincident (r ~ 20 arcsec) with several companion galax- 
ies and may represent either gas associated with the ISM in 
these galaxies or stripped galactic gas from these companions. 
The region of coolest emission (T ~ 2.5 keV) corresponds to 
the peak of the X-ray surface brightness profile and is offset 
from the optical centroid of the cD galaxy by ~ 2 arcsec. How- 
ever, this scale is comparable to the map resolution in the core 
and so may be consistent with a zero offset. As we discuss in 
Paper II, the sites of coolest X-ray emission in Abell 1068 cor- 
respond to over 90% of the emergent optical blue light from the 
starburst. We defer a detailed discussion of the implications for 
the star formation history in Abell 1068 to that work (McNa- 
mara, Wise, & Murray 2004). 

In contrast, to the northwest, the temperature rises relatively 
steeply as one moves away from the central cD increasing from 
2.5 keV to 5 keV within ~ 10 arcsec (~ 25 kpc). This rapid rise 
in temperature corresponds to a sharp decline in the X-ray sur- 
face brightness at this position as seen in Figure 1 . This position 
is also coincident with an inflection point in the azimuthally av- 
eraged surface brightness profile discussed in §4. Such temper- 
ature jumps and surface brightness discontinuities are charac- 
teristic of the "cold fronts" observed in many other clusters and 
may represent evidence for previous merger activity in Abell 
1068 (Markevitch et al. 2000; Vikhlinin, Markevitch, & Murray 
2001). We note however that the surface brightness inflection 
observed in Abell 1068 is much weaker than that typically seen 
in clusters featuring cold fronts such as Abell 2142 (Markevitch 
et al. 2000). We defer a more detailed analysis of this feature 
to a later work. 

3.5. Abundance Map 

The morphology of the abundance map in Abell 1068 is even 

more striking. Figure 9 shows the metallicity as determined 

// // 

from the MEKAL fits for the same 80 x 80 central region. 
An optical image of the R-band light (F606W filter) in Abell 
1068 was obtained from the HST archive and is overlaid for 
comparison. There is a clear pattern of enrichment extending 
along the semi-major axis of the cluster and the central cD 
galaxy. Several filamentary extensions of roughly solar metal- 
licity gas are also apparent and in all cases extend from the 
center to the position of the various companion galaxies in the 
cluster. Since the map abundances are measured in projection, 
local abundances may be somewhat higher. 

Several researchers have found a trend between increasing 
metal abundance and the presence of a cD galaxy in the central 
regions of clusters (Fukazawa et al. 2000, Irwin & Bregman 
2001, David et al. 2002). This trend may be somewhat stronger 
in cooling flow clusters (Irwin & Bregman 2001). Two of these 
studies were carried out using the low resolution ASCA and 
BeppoSAX observatories, and were based on radially-averaged 
trends between metallicity and radius. These trends suggest 
that the correlation between the presence of the cD and the cen- 
tral metallicity enhancement came about primarily as a result 



of type one supernovae (SNe la) in the cD's stellar populations, 
with some contribution of type two supernovae (SNe II) from 
star formation in the cooling flows. 

In Figure 9, we show that there is a remarkable correlation 
between the location and shape of the Abell 1068 cD galaxy's 
optical isophotes and the isometallicity contours. The strongly 
flattened appearance of both the isophotes and the metallicity 
contours renders a chance coincidence unlikely. It suggests 
strongly that the presence of the cD and regions of enhanced 
metallicity are strongly coupled. The average cluster metal- 
licity beyond the cD is Zq ~ 0.4, in solar units, and rises to 
Z c d ~ 0.8 at the location of the cD. This metallicity increase has 
likely occurred either by the external introduction of enriched 
gas during a group merger or by an internal process related to 
mass loss from the stars in the cD or the cluster at large. Given 
the correlation between the cD's optical isophotes and the con- 
tours of highest metallicity, we will examine whether mass loss 
from the general population of stars in the cD and the ongo- 
ing star formation would be capable of producing the observed 
metallicity enhancement. 

We can characterize the metallicity of the keV gas in the 
vicinity of the cD as Z cD = Zq + SZ, where SZ represents the 
metal enhancement due to supernova explosions, then 

V M H J \M h Jq v ' 

In order to calculate the supernova type la rate, we adopt Tam- 
mann's (1974) value 4 x 10~ 13 Lb yr -1 . The integrated /?-band 
absolute magnitude of the cD out to a 170 kpc radius is -24.0, 
which corresponds to L(R) = 2 x 10 n L Q . Adopting an average 
color for a typical cD of (B-R) = 1.6 (Peletier et al. 1990), we 
find an integrated B-band luminosity of L(B) = 1.4 x 10 n L Q . 
This corresponds to a SNe la rate of 0.03 yr" 1 . For the the iron 
yield per SNe la we adopt 0.74 M /SNe I (cf. David et al. 
2001), giving a total yield of y 7 - 2 x 10 7 M Q (f/Gyr). 

In addition to the yields from SNe la, there will be a sig- 
nificant contribution of metals from SNe II associated with 
the ongoing star formation in the cD. Assuming an iron yield 
of 0.12 M /SNeII, and a supernova rate of 1 per 100 M Q 
of star formation (cf. David et al. 2001), we obtain a to- 
tal yield of y n ~ 1.2x 10" 3 M A p M Q . Adopting the luminos- 
ity masses for the accretion populations in Table 2, we find 
y n ~ 8 x 10 5 - 4 x 10 6 M Q . In order to match the observational 
constraints, the supernova yields must approximately double 
the ambient metallicity of ~ 0.4 solar by seeding the ambient 
gas with metals. >From Figure 13, we find an integrated gas 
mass within a 170 kpc radius of 3.4 x 10 12 M Q . Further adopt- 
ing a solar iron abundance (M Z /M H )Q = 4.67 x 10~ 5 , we find 
that the observed central abundance gradient would accrue af- 
ter only ~ 3 Gyr. This conclusion is largely independent of the 
SNe II yields from the starburst. Assuming SNe II yields from 
the maximum mass of the accretion population permitted by the 
data, Map = 3.1 x 10 9 M (cf. Table 2), the mass ratio of SNe 
la to SNe II yields is 15. The implications of this result for the 
starburst in Abell 1068 is discussed in more detail in Paper II. 

To summarize, the metallicity gradient associated with the 
cD galaxy can be almost entirely accounted for by type la su- 
pernovae originating in the ambient stellar population of the cD 
over a reasonable age for the galaxy. However, in the central 
starburst region, the SNe II yields from the starburst will be- 
come increasingly important. Furthermore, the central metal- 
licity will almost certainly continue to grow as the roughly 
4 x 10 10 M Q of remaining molecular fuel (Edge 2001) is con- 
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sumed by star formation. We should emphasize the factors of 
several uncertainty in both the star formation masses, supernova 
rates, and yields. However, the numbers taken at face value, 
coupled with the observed correlation between the optical and 
metallicity isocontours, present persuasive case for cD galaxy- 
driven metallicity gradients in the keV gas. 

4. SURFACE BRIGHTNESS PROFILES 

Despite the apparent spatial complexity evident in the core 
of Abell 1068, the overall surface brightness in the cluster can 
still be well described by an analytic (3 model profile. This 
profile combined with the temperature information from §3.2 
can be used to determine the density profile, and various other 
physical properties of the ICM. Surface brightness profiles were 
calculated by constructing image mosaics of all chips in the ob- 
servation in three bands: a total band image from 0.3-7.0 keV, 
and soft and hard bands spanning 0.3-2.0 keV and 2.0-7.0 keV, 
respectively. Point sources were filtered out as discussed in §2. 
Appropriately weighted exposure maps were then created for 
each band in order to accurately correct for variations in the de- 
tector efficiency over the field. Finally, the surface brightness 
profiles were calculated in elliptical annuli, 1 arcsec in width 
with a position angle and axis ratio as discussed in §2. The 
profiles were determined out to a maximum radius of 16.6 ar- 
cminutes or 2.8 Mpc. 

Figure 10 shows the resulting surface brightness profile in 
the 0.3-7.0 keV band. Profiles for the soft and hard bands look 
similar. An inflection in the profile is apparent at a radius of 
~ 10 arcsec (<~ 25 kpc). Attempts to fit the surface brightness 
profiles with single (3 models yielded unacceptable fits. Conse- 
quently, we have used a double (3 model of the form 

/(r)=/ fl +/i^l + ^J + + W 

to describe the surface brightness profiles. Each component has 
fit parameters (/,-, r,-, /?,•), where the normalizations and core radii 
are denoted by /,■ and r,, respectively. The term Ig is a constant 
representing the contribution of the background. A detailed dis- 
cussion of the double (3 model and its advantages for represent- 
ing cluster X-ray surface brightness profiles can be found in 
Xue & Wu (2000) and Ettori (2000). 

Table 2 shows the results from fitting equation 2 to the pro- 
files for the total, soft, and hard bands. In all three bands, 
a broad component corresponding to the outer cluster gas is 
present with fit values of r c <~ 24 arcsecs for the core radius 
and (3 <~ 0.59, consistent within the measurement errors. For 
the total 0.3-7.0 keV, the fit to the inner component yielded val- 
ues of r c ~ 8 arcsecs and (3 ~ 0.99. As Figure 1 1 illustrates, 
these two components are distinct at the 99% confidence level. 
Comparing the soft and hard band fits, we find that the inner (3 
component in the soft, 0.3-2.0 keV band is more extended rel- 
ative to the hard band with fit parameters of r c <~ 1 1 arcsec and 
(3 ~ 1.0. This extended soft core is consistent with the flatten- 
ing of the temperature profile inside 10 arcsec discussed in §3.2 
and shown in Figure 3. >From 2.0-7.0 keV, the best fit value of 
the core radius for the inner component is r c <~ 4 arcsecs, possi- 
bly due to the presence of non-thermal emission from the unre- 
solved, central AGN or heating from the starburst. Reduced x 2 
values for the total, soft, and hard band fits were 1.8, 1.2, and 
1.2, respectively. 

We note that recent models of bubbles in cluster atmospheres 
associated with central AGN can produce "kinks" or inflections 



in the azimuthally averaged surface brightness similar to that 
seen in Figure 10. However, no other obvious evidence for 
bubbles is apparent in the images of Abell 1068. Also, the 
strong correlation between the optical isophotes and the abun- 
dance map presented in Figure 9 is difficult to reconcile with 
significant mixing of the ICM due to bubbles. 

5. PHYSICAL CONDITIONS IN A1068 

5.1. Density Profile 

Following the derivation of Xue & Wu (2000), we assume 
that the two (3 model components in equation 2 correspond to 
two phases of gas in the ICM with potentially different electron 
temperatures T\ and Ti. With this assumption, equation 2 can 
be inverted to yield the electron number densities for the two 
components as well as the combined electron number density, 
n e (r), 



n e (r) = y^n ei (r) 

i 

n ei {f) = 



1 1/2 



i 

n e (oy 



n e (r) 



n ei (r) 



n ei (r) = n ei (0) 1 + 



-3A 



, ' = 1,2 (3) 
(4) 
(5) 



where the core radii, r,, and exponents, /?,, are the same as be- 
fore and n e (0) represents the central, total electron density. The 
central number densities for the individual components are re- 
lated to the fitted surface brightness parameters, (/,-, r,-, /?,■), by 



4(0) = 



47T 1 / 2 



IX3A-) 



1X3 A- 1/2) 



where 



£-(*)(*)(*) 



1/2 



roft) r(3/3 r 



1/2) 

T72) 



(6) 



(7) 



7 = 1,2 and jji 

and gi is the Gaunt factor for component i and a(T) is the emis- 
sivity due to thermal bremsstrahlung as given in Xue & Wu 
(2000). 

The choice of values for T\ and T2 is somewhat uncertain; 
however, as equations 6 and 7 indicate, the temperature depen- 
dence of the central electron densities is fairly weak. We have 
chosen the minimum and maximum temperature values from 
the deprojected temperature profile discussed in §3.3 for the 
narrow and extended components, respectively. Using T\ = 1 .76 
keV and T2 = 5.0 keV, we obtain central densities of n e i(0) = 
0.135 cm" 3 and n e 2(0) = 0.024 cm -3 . The resulting density pro- 
files are show in Figure 12. We note that using the two temper- 
atures derived from the two-phase model (2T) fits in Table 1 re- 
sulted in central densities within 5% of these values. Of course, 
if the gas is truly multiphase as the spectral analysis in §3.1 in- 
dicates, then representing the temperature structure using only 
two components is simplistic. However, the observed surface 
brightness profile does not allow additional components to be 
unambiguously identified. 

Using the derived density profiles, the spectral fitting results 
of Table 1, and the assumption of hydrostatic equilibrium, we 
have computed various estimates for the gas and total mass in 
Abell 1068. Table 3 gives the gas mass and total mass within 
0.5 Mpc for the various spectral models discussed in §3.1 fol- 
lowing the method discussed in Hicks et al. (2002). Figure 13 
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shows the gas mass profile for the best fit FC model. In addi- 
tion to the total cluster estimates, we have calculated the mass 
associated with the central, condensed /3 model component sep- 
arately. These estimates are labeled "Core" in Table 3 and as- 
sume a central temperature of 1.76 keV taken from the depro- 
jected temperature profile discussed in §3.3. 

5.2. Cooling Time 

The characteristic time that it takes a plasma to cool isobari- 
cally through an increment of temperature ST can be written 

- £ kST 

tcOOl — r. . (8) 

2 n e A(T) 

where n e is the electron density, A(T) is the total emissivity of 
the plasma (the cooling function), and k is Boltzmann's con- 
stant. For isochoric cooling, 5/2 is replaced by 3/2. Utilizing 
the temperature profile from §3.2 and the density profile from 
§5.1, we can calculate the cooling time as a function of radius 
in the cluster. Figure 12 shows the resulting cooling time as a 
function of radius in Abell 1068. The central cooling time in 
the core reaches a minimum value of 9 x 10 7 yrs. The "cooling 
radius" where the cooling time exceeds the age of the cluster, 
canonically taken to be 10 10 yrs, is 282 kpc. As we show in 
Paper II, 98% of the blue light associated with star formation in 
Abell 1068 occurs inside the central <~ 40 kpc where the cooling 
time is < 5 x 10 8 yrs. 

5.3. Entropy 

Similarly, we can use the density and temperature informa- 
tion to calculate the entropy in the cluster gas as a function of 
position using 

S=T/n 2/3 (9) 

where T is given in keV. Figure 12 shows the resulting entropy 
profile. This figure shows what is becoming an ubiquitous fea- 
ture of cluster cores, an "entropy floor" in the central regions. 
Inside of ~ 20 kpc, the gas is essentially isentropic. Various 
mechanisms have been proposed to explain the presence of this 
entropy floor including pre-heating of the gas (Kaiser 1991) by 
supernovae-driven winds (e.g. Ponman, Cannon, & Navarro 
1999) and active galactic nuclei (e.g. Yamada & Fujita 2001; 
Nath & Roychowdhury 2002). Alternatively, it has been argued 
that radiative cooling can be an efficient means for introducing 
an entropy floor at the observed levels. Voit et al. (2002) have 
studied the effect various entropy modifications on the proper- 
ties of the cluster gas. 

6. CONDUCTION 

XMM RGS and Chandra HETG observations of cluster cores 
have failed to detect emission lines from gas at temperatures 
less than about 3 keV indicating that cooling may be truncated 
somehow (Peterson et al. 2001, Peterson et al. 2003, Wise et 
al. 2004). These data have led to the revitalization of various 
heating mechanisms previously considered to be ineffective in 
clusters including heating by the central AGN and thermal con- 
duction. Recent work by Narayan & Medvedev (2001) and Za- 
kamska & Narayan (2003) however suggests that thermal con- 
duction may be effective in clusters. Using the derived density 
and temperature in the cluster, we can directly calculate the heat 
input due to thermal conduction in Abell 1068. Following the 
analysis of Zakamska & Narayan (2003), we can write 

F = -k^- (10) 
dr 



for the heat flux due to thermal conduction, where the thermal 
conductivity is 

/1.84x 10" 5 r 5 / 2 \ _, _, _, 
K = fK sp =f\ — J ergss K cm (11) 



■ 



and In A ~ 37 is the usual Coulomb logarithm. Here, the con- 
ductivity is expressed relative to the canonical Spitzer value, n sp 
by a fixed fraction /. Figure 14 shows the resulting conductive 
heat flux as a function of radius in the cluster assuming / = 1, 
i.e. assuming thermal conduction is operating at the Spitzer 
value. The heat flux peaks at a radius of ~ 50 kpc where the 
temperature gradient is steepest. 

A more interesting quantity to calculate, however, is the re- 
quired level of thermal conduction necessary to balance cool- 
ing. If we assume that the observed X-ray luminosity due to 
cooling within a given radius is completely balanced by heating 
due to thermal conduction from hotter gas outside that radius, 
we can write 



L x «r) = Anr L f k Sj 



~d7J' 



(12) 



Solving this equation for the required value of the conduction 
parameter, /, at each position in the cluster yields Figure 15. 
As this figure shows, the thermal conductivity of the gas must 
exceed the Spitzer value by an order of magnitude to balance 
radiative cooling over the cooling radius in Abell 1068. Conse- 
quently, it seems unlikely that thermal conduction could quench 
cooling in Abell 1068, especially in the core. Similar results 
have been found for analyses of other cooling flow clusters ob- 
served by Chandra including Abell 2052, Abell 2597, Hydra A, 
Ser 159-03, and 3C 295 (Zakamska & Narayan 2003) as well 
as the cores of Abell 2199 and Abell 1835 (Voigt et al. 2002; 
Fabian, Voigt, & Morris 2002). As we show in Paper II, energy 
input from other sources such as the central AGN and super- 
novae can account for < 25% of the X-ray luminosity from 
cooling gas. 

7. CONCLUSIONS 

We have performed a detailed spatial and spectral analysis of 
a 26.8 ksec Chandra X-ray observation of the Abell 1068 clus- 
ter of galaxies. This cluster is exceptional due to the presence 
of a massive 20 — 70 M Q yr _1 starburst in the core. Our primary 
goal was to determine whether or not the data for this object are 
consistent with the standard cooling flow model where cooling 
X-ray plasma accumulates in some cold form such as stars. Our 
analysis indicates that Abell 1068 exhibits many of the common 
characteristics seen in other clusters with cool cores including 
a sharply peaked surface brightness profile, declining temper- 
ature gradient, and flat entropy profile in the core. Although 
discriminating between various spectral models can be difficult 
at CCD resolutions, the integrated X-ray spectrum for Abell 
1068 is best fit by a cooling flow model with a total mass de- 
position rate of <~ 150 M Q yr _1 . Inside a radius of 40 kpc, the 
measured cooling rate drops to ~ 40 M Q yr _1 which is com- 
pletely consistent with the star formation rates. 

Although measuring mass depositions rates from X-ray data 
is notoriously difficult to interpret due to the degeneracy in 
spectral fits, the cooling time in the gas is a well constrained 
quantity. The X-ray surface brightness distribution for Abell 
1068 provides an accurate measure of the density in the gas. 
When combined with the measured temperature profile, we find 
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that the cooling time of the gas in the core of Abell 1068 is very 
short, reaching a minimum of 9 x 10 7 yrs. The cooling time 
inside 40 kpc, where 98% of the star formation is observed in 
Abell 1068 occurs (Paper II), is < 5 x 10 8 yrs. The close spa- 
tial correspondence between regions of observed star formation 
and short cooling times again points to a connection between 
the cooling X-ray plasma and the star forming material. 

The abundance profile in Abell 1068 is especially striking. 
The annular analysis clearly indicates a central enhancement 
in the core roughly a factor of 2 higher than in the outer re- 
gions. We have also constructed a 2D map of the abundance 
in the cluster which shows a strong spatial correlation with the 
optical light. Our analysis indicates that this central increase 
is consistent with enhancement by supernova explosions asso- 
ciated with the cD's ambient stellar population and the recent 
starburst. This correlation also argues against the presence of 
significant energy input due to "bubbles" generated by a central 
AGN, since such a mechanism would tend to smooth out such 
enhancements on a relatively short timescale. And in any event, 
the central radio source in Abell 1068 is an order of magnitude 
weaker than those in clusters like Hydra A and Perseus where 
such bubbles are observed. 



Finally, we have examined the possibility that heating by 
thermal conduction could quench cooling in this object. We 
find that conduction is an order of magnitude too weak to bal- 
ance the emission from cooling over the entire cooling region. 
In Paper II, we have examined various other potential feedback 
mechanisms including the central AGN and energy input from 
supernovae. We estimate that these mechanisms could account 
for <; 25% of the cooling luminosity. Taken all together, these 
results are consistent with a picture where the extreme starburst 
properties of Abell 1068 are being fueled by a cooling ICM. 
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n jtqo+0.006 

/_ 005 
a 574+0.035 
u "> '^-0.015 


1 994. Q+156.72 
^ -18608 
1 5+103.43 

193.(t»;« 


10.96!?il 
4.38!°;^ 


990+°- 198 

1 A65+0.333 
1 .ODJ_ 278 

554+o.in 

^•JJt_0.O55 


77 1 7 4+527.77 
" 1 ' -^-482.84 

5950.4^4.37 
1377.6^ 


7 547+0.04 
Z.Jtz,_oo4 

1.050® 
1.458t^ 


1.77 
1.19 
1.19 


723.7 
486.0 
488.6 


417 
416 
417 



a Surface brightness / in units of 10 9 photons sec 1 cm 2 arcsec 
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Table 3 
Mass Estimates 


Model 


kT 
[keV] 


no 

[10- 2 cm- 3 ] 


M gas 
[1O 13 M ] 


[10 13 M Q ] 


IT 
2T 
FC 
TC 
Core 


o 76+0.22 
J - /o -0.14 

6.85^ 

4 91+0.21 

£T Q9+O.66 
J - yz -1.01 

i - /o -0.10 


9 40 +o.i7 

z - w -0.17 
9 QO+0.17 

9 44+0.17 
Z -^M).18 
9 07+0.17 
^""-0.17 

13.53!j;| 


4 79+0.90 
H ' -0.80 

4.68S 

^•°'-0.79 

1 s +0 - 32 

u - i;5 -0.10 


25.28!}| 

4c 4^+12.08 

27.41+f 3 ' 
41.21+™ 3 
19.77^ 
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FIG. 1 . — An exposure corrected, 0.3-7.0 keV image of the central 80 X 80 of Abell 1068 from the 26.8 ksec ACIS-S exposure. The image has been adaptively 
smoothed using the CIAO tool CSMOOTH. The optical contours corresponding to the HST R-band (F606W) surface brightness are also shown. Note the offset 
between the X-ray centroid and the optical center of the cD galaxy. 
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FIG. 2. — Integrated 0.3-7.0 keV spectrum for Abell 1068 from the 175x129 arcsec elliptical extraction region discussed in the text. The spectrum has been 
adaptively binned to contain 30 counts per bin. The solid line shows the best fit spectral model, in this case the FC cooling flow model discussed in the text. 
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FIG. 3. — The projected gas temperature for Abell 1068 as a function of projected radius along the semi-major axis in the cluster. One sigma errors bars are 
indicated. The horizontal dashed line marks the mean gas temperature from the best fit to the integrated spectrum shown in Figure 2, while the hatched horizontal 
region depicts the 90% confidence interval for the temperature. The solid line shows the best fit analytic model. 
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FIG. 4. — The projected abundance profile for Abell 1068 as a function of projected radius along the semi-major axis in the cluster. One sigma errors bars are 
indicated. The horizontal dashed line marks the mean abundance from the best fit to the integrated spectrum, while the hatched horizontal region depicts the 90% 
confidence interval for the mean abundance. 
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FIG. 5. — The best fit foreground absorbing column in Abell 1068 as a function of projected radius along the semi-major axis in the cluster. One sigma errors 
bars are indicated. The horizontal dashed line marks the column from the best fit to the integrated spectrum, while the hatched horizontal region depicts the 90% 
confidence interval for the mean column. 
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FIG. 6. — The enclosed mass deposition rate for Abell 1068 as a function of projected radius along the semi-major axis in the cluster. One sigma errors bars are 
indicated. The horizontal dashed line indicates the total deposition rate from the best-fit FC cooling flow model discussed in the text. The vertical dotted line marks 
the cooling radius as derived from the temperature and density profiles. Two analytic forms for Mx(< r) are depicted by the dot-dashed and dot-dot-dot-dashed 
lines: the canonical Mx(< r) oc r, which does not describe the profile well, and Mx(< r) oc r 0A , respectively. 
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FIG. 7. — The deprojected gas temperature for Abell 1068 as a function of projected radius along the semi-major axis in the cluster. One sigma errors bars are 
indicated. The horizontal dashed line marks the mean gas temperature from the best fit to the integrated spectrum shown in Figure 2. 
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Abell 1068 




FIG. 8. — A map of the gas temperature for the central 80 X 80 in Abell 1068. The overlaid contours show the surface brightness from the adaptively smoothed 
flux image shown in Figure 1. The color bar on the right indicates the temperature scale in keV. Errors for the measured temperatures range from 10% — 30% over 
the map. 
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FIG. 9. — A map of the elemental abundance for the central 80 X 80 in Abell 1068. The optical contours corresponding to the HST R-band (F606W) surface 
brightness are shown. The color bar on the right indicates the abundance in units of solar abundance. Errors for the measured abundances range from 10% — 30% 
over the map. 



Chandra Observations of Abell 1068 



21 



a1068_dl_290_7000_filt_fit.ps 




R [arcsec] 

FIG. 10. — Radial surface brightness for Abell 1068 in the 0.3-7.0 keV band accumulated in 1 arcsec annular bins. The solid lines indicate the two components 
of the best fit double (3 model discussed in the text. The horizontal dotted line represents the best fitting background value. The lower pane shows the residual 
deviations from the fit. 
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FIG. 1 1 . — Confidence contours for the two components of the double j3 model fit to the surface brightness profile shown in Figure 10. Contours are shown for 
la, 2a, and 3a. The crosses indicate the location of the best fit (r c , (3) values for the two components. 
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FIG. 12. — Radial variation of (a) electron density, (b) isobaric cooling time, (c) pressure, and (d) entropy. In panel (a), the dashed lines depict the two components 
of the electron density from the double /3 model fit to the surface brightness distribution. In panel (b), the horizontal line indicates a cooling time equal to 10 years 
and the vertical line indicates the corresponding cooling radius. 



24 



Wise, McNamara, & Murray 




Fig. 13.— 



Gas mass profile for Abell 1068. 
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FIG. 14. — Heat flux due to thermal conduction as a function radius in the cluster assuming conduction operates at the Spitzer value (i.e., k = K sp ). 
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FIG . 15. — Value of the conduction parameter, f, required to balance radiative cooling as a function of radius in the cluster. The conduction parameter is expressed 
as a ratio relative to the Spitzer value. 



